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ABSTRACT 

We estimate effects of the final state interactions in B — > tttt decays coming from rescattering 
of tttt via exchange of p, a, /o mesons. Then we include the pp rescattering via exchange of tt, lj, 
ai mesons and finally we consider contributions of the aiir rescattering via exchange of p. The 
absorptive parts of amplitudes for these processes are determined. In the case of 7r + 7r~ decay 
mode, due to model uncertainties, the calculated contribution is \Ma\ < 1-7 x 1CT 8 GeV. This 
produces a small relative strong phase for the tree and color-suppressed B — > tttt amplitudes 
consistent with the result of a recent phenomenological analysis based on the BaBar and Belle 
results for the B — > tttt branching ratios and CP asymmetries. 



1 INTRODUCTION 



The experimental results on B decays coming from Belle and BaBar offer many puzzles for 
theoretical studies. Among them the B — > tttt decays are particularly interesting Q]I2|- Many 
theoretical frameworks such as perturbative QCD approach of Beneke, Buchalla, Neubert and 
Sachrajda (BBNS) |3I and the approach of 4 , Soft Collinear Effective Theory (SCET) [5- 
9] and many others [10-24] have attempted to understand the observed decay rates. Within 
QCD factorization charmless two body decays of B mesons have amplitudes which factorize at 
lowest order in 1/mt,. It means that in this approach, in neglecting the next-to- leading terms 
in l/m& expansion, one ends up with the naive factorization ansatz. The naive factorization 
(e.g. jlUl lllj ^ gave the rate of B° — > 7r + 7r~ too large in comparison with the observed rate 
while the B° — > 7r°ir° decay rate came out too small within this simple framework. Agreement 
with experimental data on5-> tttt has been found within both BBNS and SCET frameworks. 
The improved B — ► vr 7r° decay rate was obtained recently within BBNS |Hj with the presence 
of parameter A& whose precise value is unknown 21 . Within SCET the agreement with the 
experimental data is achieved [2] with the presence of non negligible long-distance charming 
penguin contributions. It has been pointed out in Ref. |25j that in B weak decays one cannot 
neglect the effects of final state interactions due to the growth of forward scattering of the final 
state with the squared center off mass energy, as required by the optical theorem and cross 
section data. This indicates that "soft scattering does not decrease for large tub" [2*5] . 

Recently the authors of [21] considered two-body decay modes by including final state inter- 
actions (FSI). Contributions of the cc state, which in the literature very often called charming 
penguins were considered in EE]- The charm meson rescattering due to charm meson ex- 
change has been considered in Refs. |131 I23| and more recently in [21J. It was found the largest 
contribution appears in the B — » Kit mode |13| , but is much smaller in the case of tttt final state 
|23| . The authors of k 24 found that the absorptive part of the rescattering cannot explain the 
observed enhancement of the branching ratio and cannot produce a small branching ratio 
of the 7r + 7r~ rate. 



framework as described in [21], but we point out that there are more intermediate states which 
contribute to both amplitudes and give important contributions. As in [21] we take into account 
only dominant contributions proportional to the effective Wilson coefficient a\. In this approach 
for the charmless final state interactions only the contributions of tttt and pp intermediate states 
were used in [2%] . Since in B decays, resonant FSI is expected to be suppressed due to the 
absence of resonances at energies close to the mass of the B meson, we consider only t- channel 
FSI. However, in the case of tttt — ► tttt rescattering we include possibility that in addition to 
the p meson exchange there are contributions coming from a and /o exchange. In the case of 
pp — * tttt rescattering we find that there is a contribution of the u meson for the tt + tt~ final state 
as well as contributions of the ai(1260) axial meson. We determine contributions coming from 
ai(1260)7r intermediate states, inspired by the recent BaBar measurement of the very large rate 
for B° -> aj~7r+ state with the branching ratio BR(B° -> o|(1260)vr-) = (40.2±3.9±3.9) x 1CT 6 
|27| . In our approach the ai(1260)~7r + rescatter via p° exchange into the 7r + 7r~ final state. 
Although the B° — > afir~ decay rate has not been observed yet, we estimate this contribution 
assuming the naive factorization for the amplitude. The paper is organized as follows: in Sec. 2 
we give basic formulas for the two-body B amplitudes and the Lagrangian describing the strong 
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interactions of the light mesons used in our calculations, in Sec. 3 we present results of our 
calculations for the absorptive part of the amplitude, in Sec. 4 we discuss our results and we 
summarize them in Sec. 5. 



2 THE FRAMEWORK 

In the studies of the B — > tttt and B — > Kir branching ratios and CP asymmetries it was 
found that amplitudes arise from tree, color-suppressed, penguin and the electroweak penguin 
diagrams (see e.g. |191 124) ). In our approach we consider only leading contributions in charmless 
FSI and therefore we only use the effective weak Lagrangian for the process b — > udu at the tree 
level in the following form: 



G 

C w = --j=V ub V* d ai{ub) V - A {du) v ^ A - 



Here a\ is the Wilson coefficient and we use the same value as given in |24j (ai(/i) = 0.991 + 
i0.0369; the scale fi = 2.1 GeV), which includes short-distance nonfactorizable corrections such 
as vertex corrections and the hard spectator interactions determined within QCD factorization 
approach [3] . In our further study we use naive factorization approximation jlUj , in which the B 
meson decay amplitude can be written as a product of two weak current matrix elements. The 
standard decomposition of the weak current matrix elements is: 



(V(k,e,m v )\qT»q\P(p,M)) 
e-q 



2V(q 2 ) 
M + my 



+ 2im v ^-q' 1 A (q 2 



i{M + m v ) 



M<f) 



e-q 



M + my 



M 2 



m 



A 2 (q 2 ) - 



(2) 



Similarly, heavy pseudoscalar to light pseudoscalar transition is described by the matrix element: 



(P(k,mp)\qT^q\P{p,M)) 



«2 q 



F + (q 2 ) + 



(M 2 



m. 



-«(g 2 



(3) 



while for the heavy pseudoscalar to light axial vector transition, we use the expression given in 



{A(k, e, m A )\qWq\P{p, M)) = i[(M + m A )e"V 1 (q 2 



e ■ 



2m- 



e ■ 



'Y(V 3 (q 2 )-Vo(q 2 ))]-e^e„p a } 



M + m A 
2A(q 2 ) 



P^q 2 )- 



(4) 

q- M + m A ' w 

with V 3 (q 2 ) = (M+M a )/(2m A )Vi(q 2 )-(M-m A )/(2m A )V 2 (q 2 ). In above equations <^ =p^-k^ 
and P^ = p^ + k^. The light meson creation (annihilation) is described by the matrix elements: 

(P(p)\q^(l - l5 )q\0) = ifpp» , (V(j>,e)\q 7 »(l - 75 )<?|0) = f v m v e» , 
(A(p,e)\q^(l - l5 )q\0) = f A m A e» . (5) 

In our numerical calculations we use the following values of relevant parameters as given in [21] : 
U = 0-132 GeV, f p = 0.21 GeV, f ai = 0.205 GeV, F B *(0) ~ F^{ml 



0.25 ~ F^(m 2 ai ), 
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A Bp (0) ~ A Bp {m 2 p ) = 0.27, A Bp {0) ~ A Bp (m 2 p ) = 0.26. We use: V Bal (0) ~ ^ Bol (ro*) = 0.13 

ESI- ' 

Using above expressions the leading contribution to the amplitude for B° — ► 7r 7r + was found 
to be (e.g. dU) 

^(5° - tt+O = = -i^V^V: d ai[F^(ml)(ml - ml)]U . (6) 

In j2U the value ai = 0.9921 + i0.036 led to the amplitude A*{B -► 7r + 7r-) SD = 3.2 x 10~ 8 + 
il.2 x 10~ 9 GeV (we took the V u b = 0.00439 [25])- Without color-suppressed and penguin 
contributions this gives the branching ratio BR(B° — > vt + 7t~)sd = 9 x 10~ 6 , too large in 
comparison with the average experimental value (4.6 ±0.4) x 10~ 6 as given in [21]. The inclusion 
of color-suppressed and penguin amplitudes decreases the rate ^2I2I]> but it is still too large 
in comparison with experimental result. 
The amplitude for B° — ► p + p~ is 

/ —2iV(m 2 ) 
A(B°(pb) - P + (qi, e 1 )p-(q 2 , e 2 )) = U p e^e^q^gsp 



M B + m p 



+A 1 (m%M B + m p ) £l • e 2 - 2A 2 (m 2 p ) ei - PbPb ' " 2 \ , 

Mb + m p J 

with A p = —-^VubV* d aif p m p . The amplitudes for B° — ► a^7r + and S° — > a^7r~ are: 



(7) 



A(B°(p) -» O! fe, e)vr + (gi)) = iA ai ,i(p + qi) ■ e, 

A(S°(p) -» a+(gi,e)vr-(g 2 )) = i4 0l , 2 (p + gi) • e, (8) 

with A ai ,i = -^V ub V: dai f ai m ai F^(m 2 al ) and Au,2 = -^V^V* d a x U2m ai V Q Ba \ml). 
The light mesons' strong interactions are described by 

Strong = i^=-Tr(p^[U, d p U}) - 4^j^e p ^Tr(d p p u d a p p Ii) 

+ G A vpTr(A p [p»,n]) + iG s V2Tr(nUS) + iG s >V2Tr(ILILS') . (9) 

In these equations II is the 3x3 matrix containing pseudoscalar mesons, p is the 3x3 matrix 
describing light vector mesons, and S, S' are matrices describing scalar mesons. In our numerical 
calculations we use g pnn = 5.9 and Cyyp = 0.33 (see [SD] - |32]). The coupling |GvLKp| = 3.12 
GeV is obtained from the experimental results for a® — ► p~n + decay width Ta = 0.2 GeV. 
Finally, the couplings G s and G' s are obtained by using PDG data 1 on a (or /o(600)) and 
/ (980) meson: m a m (0.4-1.2) GeV, T a m (0.6-1) GeV, m f = 0.98 GeV and T j « (0.04-0.1) 
GeV. In the numerical calculation we take the average values m a = 0.8, T(a — ► 7T7t) = 0.8 GeV, 
m/ = 0.98 GeV and r(/ (980) tttt) = 0.07 GeV and we determine G s = 4.24 GeV, and 
G' s = 1.37 GeV. 

Using naive factorization we obtain for the branching ratio BR(i?° — > a^7r + ) = 1.8 x 10~ 5 
about two times smaller than the experimental result given in |27| . Using above mentioned data 
we predict that BR(£° -> a^TT - ) = 8.2 x 10~ 6 . 
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Figure 1: Feynman diagrams for Bq — > n + n decay coming from rescattering of tttt via exchanges 
of p, a, fo , pp rescattering via exchanges of w, u, a\ and a\ir rescattering via exchange of p. 

3 THE ABSORPTIVE PARTS OF THE AMPLITUDES 

In our calculation of the absorptive parts of amplitudes we include the contributions coming 
from the graphs presented in Fig. 1. The absorptive parts of amplitudes are obtained when the 
cut is done over the intermediate states tttt, pp and a\Ti as schematically given in Fig. 2. In our 
further formulas we denote momenta of particles as given in Fig. 2. The couplings describing the 
strong interactions of light mesons in these diagrams are all far of mass shell. In the approach 
of |241 I33| the additional form factor was included. Its role is to take care of the off-mass shell 
effects [HH: 

= (£fi§- < 10 > 

where t{y) = {q\ — k\) 2 , A = M3 + Aqcd and M3 is the mass of exchanged particle A3 (see Fig. 
2). We take Aqcd = 0.3 ± 0.05 GeV. Following the contributions given in Fig. 1, we determine 
the absorptive parts of the amplitudes 



" pwn 327rm| J-i y wKyj 2ml - 2S(y) - mj 
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Figure 2: The absorptive parts of amplitudes are obtained when the cut is done over the inter- 
mediate states A\ and Ai- 



M A --AAG S j^ dy — (12) 
MY f ~ -A iC >2 ^ 1/2 (rn 2 B ,ml,™l) f l , F*(y,m f ) 



^ 32vrm| J _, ay m l + mj - 2S(y) - ml 

((m B + m^^m 2 )^^) - 2A 2 (m2)/(m B + m p )C p , 2 (y)) , (14) 

Ma ~ Ap { U ) 32vrm| U dV m% + mj - 2S(y) - ml 



((m B + m^^ra^C^d/) - 2A 2 (m 2 p )/{m B + m p )C pA {y)) , (15) 
M P Pai _ A C 2 X 1/2 (m 2 B ,m%ml) yi F 2 (y,m ai ) 

" 2 "VwP — -y J 1 d y ml + m 2_ 2S{y) _ m 2 i X 

((m B + m^ij^)^^!/) - 2A 2 (m2)/(m B + m p )C P)6 (y)) , (16) 



■M^ = -Aii,lV2GrAyPfl , /OT 7r 327rm 2 



01 ' X 



/>^>S3^^' (17) 



/' *c.,A)-5- rill 2. ( 18 > 

where with C stands for the functions of momenta defined in Appendix A, while S(y) is the 
scalar product: 

S{y) = h ■ gi = hoE! - (19) 
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A = 0.25 GeV 


A = 0.30 GeV 


A = 0.35 GeV 


" VP/ 
7T7r(cr) 

7T7r(/ ) 


13 3 + 5? 
-0.5 - 0.02? 
-0.03 - 0.001?' 


17 3 + 6? 
-0.6 - 0.02? 
-0.05 - 0.002?' 


21 8 + 8i 
-0.8 - 0.03? 
-0.06 - 0.002?' 


V 

^7T7r 


12.5 + 0.5i 


16.7 + 0.6? 


21 + 0.8? 


pp(vr) 

pp(m) 


-1.7-0.06? 

5.5 + 0.2? 
-0.9 - 0.03?' 


-2.2 - 0.08?' 

7.7 + 0.3? 
-l.4-0.05i 


-2.8-0.1? 
10.3 + 0.4?' 
-1.6-0.06? 




2.8 + 0.1?' 


4.3 + 0.2?' 


5.9 + 0.2?' 




5.6 + 0.2?' 
1.9 + 0.1? 


7.5 + 0.3?' 
2.5 + 0.1» 


9.5 + 0.3?' 
3.2 + 0.1?' 



Table 1: The absorptive parts of amplitudes coming from the diagrams M\ X 10 7 V^[GeV] 
given in Fig. 1. 

and y = cos(kuqi). We use |cfi| 2 = -^rX(m 2 B , Mf, Mf), |A?i| 2 = ^A(m|, m 2 , m 2 ) and 

Ei = I 1 2 + -^i an d &io = l^i 1 2 + m2 - Here iW, stands for the masses of intermediate particles 
Ai and A(a, 6, c) = a 2 + 6 2 + c 2 — 2a6 — 2cfo — 2ac as usual. 

4 DISCUSSION 

After numerical evaluation 1 of these integrals we present our results in Table 1. We give values 
of the absorptive parts of the amplitude for three different values of the scale A = 0.25; 0.3; 0.35 
GeV. As seen from the table these amplitudes are sensitive to the choice of this parameter. It is 
important to note that the relative sign of these contributions cannot be completely determined 
|33| . By assuming that strong couplings do not have any phases, the sum of contributions 
coming from the irir —* tttt rescattering is then E T7r = (1.7 + 0.06?) x 10 -6 V^, GeV, which for 
\V ub \ = 0.00439 gives |S^| = 7.5 x 10~ 9 GeV (for A = 0.3 GeV) . It is interesting that the 
exchanges of scalar mesons give very small contributions. The contribution of pp intermediate 
states with the exchanges of vr°, u and ai is about four times smaller than the total ir + ir 
intermediate state contribution. Among these the effect of the lo exchange is important. This 
contribution was not considered in [21]. The contributions of aiir intermediate states might be 
significant, close in size to the leading tt + tt~ elastic-rescattering effect. Then in the best case (by 
summing the contributions given in Table 1, all with the positive signs) we can give an upper 
value for the absorptive part of the amplitude (A = 0.3 GeV): 

\M A {B° ^7r + 7r~)| < 1.7 x 10~ 8 GeV. (20) 

This value is very close in size to the short distance amplitude discussed in [21] (Eqs. (5.14)). 
On the other hand, for the certain choice of the strong couplings phases, the calculated contri- 
butions might almost cancel each other, leading to the disappearance of the absorptive part of 
FSI amplitude. 

In the case of B° — > 7r°7r° the absorptive part of amplitude comes from the same FSI and 
the upper bound is \M A (B° 7r vr )| < 1.4 x 10" 8 GeV, (A = 0.3 GeV). Note that there are no 

1 Numerical results were obtained with the help of the computer program FeynCalc |35|. 
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contributions coming from the exchanges of neutral mesons as a, fo in the case of 7r + 7r~ — > irir 
and oo in p + p~ — > 7r°7r° mode. Comparing this result with short distance amplitude given in 
|24| (Eqs. (5.14)) we see that the effect we discuss might enhance the amplitude by a factor 
of 2. However, the corresponding branching ratio is still too small in comparison with the 
experimental result. 

In order to estimate the effects of this leading FSI contribution in B~ — » 7r~7r decay ampli- 
tudes one can rely on the isospin relation 2 

A(B° -> tt+O - A{B° 7r°7r°) = -y/2A(B~ -> 7r°0- (21) 

We find that the absorptive part from tttt (elastic rescattering) and quasi-elastic FSI pp via 
the i-channel it, a±, w-exchange contributions might be important for B — * tttt amplitudes. Here 
we point out that the absorptive part of the B° — ► 7r + 7r~ amplitude produces the phase of the 
tree amplitude of [H3 EHJ while the absorptive part of the B — ► vr vr amplitude determines the 
color-suppressed phase of the amplitude in [371 EH] . In a recent paper [37j it was shown that it 
is possible to determine the strong phase separately for the tree, color-suppressed, and penguin 
amplitudes from the current BaBar and Belle measurements on B —* tttt branching ratios and 
CP asymmetries. The results show that the relative phase between the tree and color-suppressed 
amplitudes 5t — &c is rather small. Since we found the strong phase coming from calculated 
FSI effect for it + it~ (tree amplitude) and 7r°7r° (color-suppressed amplitude) to be almost of the 
same size, we can confirm the results of the phenomenological study given in Ref. 37 . 

Our calculations contain only information on the absorptive part of amplitudes indicating 
sources of uncertainties. One can in principle determine the dispersive parts of amplitudes, but 
due to many uncertainties we do not pursue in calculating these effects. As noticed in ^1 EH 
these contributions are expected to be of similar size as the absorptive parts of amplitudes for 
both 7r + 7r~ and 7r°7r° decay modes. 

Recently the authors of Ref. [20] estimated the effects of final state interactions using 
the Regge model. This analysis shows that the long distance charming penguins do not play 
important role. However, the long distance effects due to the light meson rescattering are very 
important in obtaining correct rates for B — * tttt decays }2()j . in agreement with the result of 
our calculation. 

In Ref. |3H]) using the SU(3) symmetry relations, it was found that in B — ► tttt decays 
the ratio of the color-suppressed and tree amplitudes is very large. Our calculations, obtained 
within a very different framework, confirm this finding. 

5 SUMMARY 

We can briefly summarize our results: 

1) The absorptive parts of amplitudes in B — ► tttt decays are calculated using the rescattering 
of tttt via exchange of p, a, /o; pp rescattering via exchange of it, to, a\ and contributions of the 
aiir rescattering via exchange of p. 

2) Although our results suffer from many uncertainties due to unknown relative phases and 
the dependence on the parameter A, we can say that our study shows the importance of the 

2 Note that we have used the Feynman diagram convention for the 7r°7r° amplitude as in |36|. 



7 



charmless final state interactions in B — ► tttt decays. Both the B° — > 7r + 7r and B° — > 7r 7r° 
amplitudes might get significant contributions from absorptive parts of the FSI amplitudes. 

3) Our result shows that the relative phase between the tree and color-suppressed amplitude 
5t — 5c is rather small and in agreement with the results of previous phenomenological studies. 
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6 APPENDIX A 

The functions of momenta C are (momenta fa and q are defined in Fig. 2): 

CV(y) = (?i + k 1 ) a (q 2 + k 2 f(-g a p + q a qp/m 2 ) = 2 (m\ - m\ + S(y)) , (22) 



C pA (y) = {-g a p + qi a qip/m 2 p ){-g aS + q% q S 2 / m 2 ) (2fa - Ql f(2k 2 - q 2 ) 5 

= A (-2m^mJ + m 2 B m 4 p + 2S(y) (s(y) - m|) m 2 p + m 2 B S(y) 2 ) , (23) 

C Pt2 (y) = (fa + k 2 ) a (fa + fc 2 ) 7 (-^ + qfq^/m 2 ) 

(2fci - qiH-g lS + qU/m 2 p )(2k 2 - q 2 ) s = (mj - S(y)) 2 , (24) 

Cp,3(y) = (~flW + q2aq2cr/m 2 )(-g a(T ' + QlaOlcr 1 1™%) 

(-5 7 y + (fci - gi) 7 (fci - gi) 7 //m^)e w ^e' eVp '' y, g2/ S 9i«'(*i " Qi) P (fa ~ Ox) ft 

= ml (m 2 B - 2m 2 p ) + m 2 m 2 B + 2S(y) (s(y) - m 2 B ) , (25) 



C P ,A{y) = (fa + k 2 ) a ' \-g a 'a + q2a>q2a/m 2 p )(fa + k 2 ) a (-g a<J > + qiaqia'/m 2 ) 
(-5 7 y + (fa - qi ) y (fa - <Zi)y/m^)e KCTp7 e KVp ' 7 Wgw(fci " Qi) P (fa ~ Qi)p> 

2 

= ^((m 2 B - Am 2 p ) ml + m 2 p m\ - 2S(y) (m 2 B - 2S(y))) , (26) 

Cp, 5 (y) = {-9an + q2 a q2 1 /m 2 p )(-g lp + q p q i lm, 2 a ^)(-g a p + q lp q^/m 2 p ) 
= ~ 4mjm 2 al ( 2m P ( 6m2al + " 4m '^ 

-4m^ (m^ml + Sfo) (m| - ) + m^m* + 2m|,S(y) 2 ) , (27) 
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c P ,6(y) = ( kl + ^2)2/3(^1 + k 2 ) a {-g ai + q2 a q2~ 1 /m 2 p ){-g~ 1p + q p q< fm 2 aX ) 

2 

(-<?? + qi P q?/m 2 p ) = -g^5T { 2m t ( 4 ™ai + m%) 

-2m 2 B m 2 p (3m 2 al + 2S(y)) + m 2 al m% + 2m 2 B S(y) 2 ) . (28) 



C ai ,i(y) = ( 2 <?1 + q2) a (-ga/3 + q2 a q2p/m 2 1 )(-g l3S + qpq&lm 2 ){q x + k x ) & 



+ ™4 ( 2 - s '(y) - 3m l - 2m li) + m ii 



+2m 2 B 



m% 



S(y) 



m 2 al (3m 2 B + 2S(y) 



(29) 



C ai ,2(y) = ( 2 qi + q2) a {-g a (3 + qi a qi(3/m 2 ai ){-g P& + qpq8/m 2 p ){q 2 + k 2 )s 



ir ( m l + m l ( s (y) 



2m% - 2m z al ) + 



2<i 

+m|j - m|5(y) - (m| + S(yj) ) . (30) 
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